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Abstract. An experimental study on the interaction of heavy, highly charged, and energetic ions (52 MeV 58 Ni 13+ ) with pure 
H2O, pure CO2 and mixed I-^OiCCh astrophysical ice analogs is presented. This analysis aims to simulate the chemical and the 
physicochemical interactions induced by heavy cosmic rays inside dense and cold astrophysical environments such as molecular 
clouds or protostellar clouds. The measurements were performed at the heavy ion accelerator GANIL (Grand Accelerateur 
National dTons Lourds in Caen, France). The gas samples were deposited onto a Csl substrate at 13 K. ln-situ analysis was 
performed by a Fourier transform infrared (FTIR) spectrometer at different fluences. Radiolysis yields of the produced species 
were quantified. The dissociation cross sections of pure H 2 and CO2 ices are 1.1 and 1.9 xl0~ 13 cm 2 , respectively. In the case 
of mixed FLOiCCh (10:1) the dissociation cross sections of both species are about 1 xl0~ 13 cm 2 . The measured sputtering yield 
of pure CO2 ice is 2.2 xlO 4 molec ion -1 . After a fluence of 2-3 xlO 12 ions cirr 2 the C0 2 /CO ratio becomes roughly constant 
(~0. 1), independent of the of initial CO2/H2O ratio. A similar behavior is observed for the H2O2/H2O ratio which stabilizes at 
0.01, independent of the initial H2O column density or relative abundance. 

Key words, astrochemistry - methods: laboratory - ISM: molecules - Cosmic-rays: ISM - molecular data - molecular pro- 
cesses 



1. Introduction 

H2O and CO2 are most abundant constituents of icy grain man- 
tles in the interstellar medium (Whittet et al. 1996; Ehrenfreund 
& Charnley 2000). Following Gerakines et al. 1999, a signifi- 
cant fraction of interstellar solid CO2 exists in mixtures domi- 
nated by H2O in both quiescent cloud and protostellar regions. 
In those regions the relative abundance of CO2 with respect to 
H 2 ranges from 4% to 25%. 

Inside the solar system the presence of H2O and CO2 ices 
is also ubiquitous. For example, CO2 is widely detected in 
cometary ices. Its abundance, relative to H2O, is about 6-7% in 
the comets Hale-Bopp (Crovisier et al. 1997; Crovisier 1998) 
and Hyakutake (Bockelee-Morvan 1997). H 2 and CO2 ice 
features have been observed also in the IR spectra of the icy 
Galilean satellites Europa, Ganymede, Callisto (Carlson et al. 
1999; McCord et al. 1998), at Triton - the largest moon of 
Neptune (Quirico et al. 1999) - and at surface of Mars (Herr 
& Pimentel 1969; Larson & Fink 1972). 
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Deep inside dense molecular clouds and protostellar disks, 
as well as at the surfaces of solar system bodies surrounded by 
thick atmospheres, the frozen compounds are protected from 
stellar energetic UV photons. However, X-rays and energetic 
cosmic rays can penetrate deeper and trigger molecular disso- 
ciation, chemical reactions and evaporation processes. In the 
case of solar system ices without significant atmospheres such 
as Europa, Enceladus and Oort cloud comets, the material pro- 
tection to UV photons is supplied by the upper layers (tens nm) 
of ice and only the energetic particles reach the inner layers 
(e.g. 52 MeV Ni penetrate roughly 20 fim in pure water ice). 
These ices are directly exposed to stellar photons and comics 
ray, solar wind (~ 1 keV/amu ionized particles), energetic solar 
particles and/or planetary magnetosphere particles. 

Laboratory studies and astronomical observations indi- 
cate that photolysis and radiolysis of such ices can create 
simple molecules such as CO, CO3, O3, H2CO3 and H2O2 
and others more complex organic compounds such as formic 
acid, formaldehyde, methanol, etc. (e.g. Gerakines et al. 2000; 
Moore and Hudson 2000; Brucato et al. 1997, Wu et al. 2003). 
By comparing laboratory data (production rates, formation 
cross section and half-lives for these compounds) with astro- 
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nomical observation, we will better understand the physico- 
chemical processes occurring in the astronomical sources. 

In this work, we present infrared measurements of two dif- 
ferent mixtures of H20:C02 ices as well as of pure water and 
C0 2 ices, irradiated by 52 MeV 58 Ni 13+ . In section 2, the ex- 
perimental setup is briefly described. The infrared spectra of 
frozen sample measured as a function of ion fluences are the 
subject of section 3. In section 4 a discussion is made about 
the dissociation cross sections and other radiolysis-induced ef- 
fects, as well as on the formation of new species on ices. 
Astrophysical implication is discussed in section 5. Final re- 
marks and conclusions are given in section 6. 

2. Experimental 

To simulate the effects of heavy and highly ionized cosmic rays 
on astrophysical cold surfaces, the facilities at the heavy ion ac- 
celerator GANIL (Grand Accelerateur National dTons Lourds, 
Caen, France) have been used. The gas samples (purity superior 
to 99%) were deposited onto a Csl substrate at 13 K through 
a facing gas inlet. The sample-cryostat system can be rotated 
180° and fixed at three different positions to allow: i) gas depo- 
sition; ii) FTIR measurement; and iii) perpendicular ion irradi- 
ation. 

52 MeV 58 Ni 13+ ion projectiles impinge perpendicularly 
onto H 2 0:C0 2 ice (1:1), H 2 0:C0 2 ice (10:1), as well as onto 
pure H 2 and C0 2 ices. The target ionizing effects of Ni and 
Fe projectiles, at the same velocity, are almost identical since 
they have almost the same atomic number (Seperuelo Duarte et 
al. 2010). The incoming charge state 13+ corresponds approx- 
imately to the equilibrium charge state after several collisions 
of 52 MeV Ni atoms (independent of the initial charge state) 
with matter (e.g. Nastasi et al. 1996). 

In-situ analysis was performed by Nicolet FTIR spectrom- 
eter (Magna 550) from 4000 to 600 cirT 1 with resolution of 

1 cm" 1 . The spectra were collected at different fluences up to 

2 x 1 1 3 ions cirT 2 . The Ni ion flux was 2 x 1 9 cm" 2 s" 1 . During 
the experiment the chamber pressure was below 2 x 10" 8 mbar 
due to the cryopumping by the thermal shield. Further details 
are given elsewhere (Seperuelo Duarte et al. 2009; Pilling et al. 
2010). 

Assuming an average density for the pure water ice and 
water-rich ices samples of 1 g/cm 3 and 1.8 g/cm 3 for the pure 
C0 2 ice, the thickness and the deposition rate are determined 
by using the initial molecular column density of the samples. 
For H 2 0:C0 2 (1:1) and H 2 0:C0 2 (10:1), the sample thick- 
nesses were 0.6 and 0.4 fim. The deposition rates were roughly 
7 and 3 /im/h, respectively. In the case of pure H 2 ice and pure 
C0 2 ice, the thicknesses were about 0.5 and 0.4 fim and, the de- 
position rates were 4 and 3 //m/h, respectively. The penetration 
depth of 52 MeV Ni ions is much larger than ice thicknesses 
therefore the ions pass though the target with approximately 
the same velocity (cross sections are constant). 

The molecular column density of a sample was determined 
from the relation between optical depth r„ = ln(/o/^) and the 
band strength, A (cm molec" 1 ), of the respective sample vibra- 
tional mode (see d'Hendecourt & Allamandola 1986). In this 
expression, / and 7 is the intensity of light at a specific fre- 



quency before and after passing through a sample, respectively. 
Since the absorbance measured by the FTIR spectrometer was 
Absy = log(Io/I), the molecular column density of ice samples 
was given by 

if 23 r 

N-— Tydv = — Abs v dv [molec cm" 2 ], (1) 
A J A J 

where Abs v = ln(/ //)/ ln(10) = r v /2.3. 

The vibrational band positions and their infrared absorp- 
tion coefficients (band strengths) used in this work are given in 
Table □ 



Table 1. Infrared absorption coefficients (band strengths) used 
in the column density calculations for the observed molecules. 



Frequency 


Wavelength 


Assignment 


Band strength (A) 


Ref. 


(cm" 1 ) 


(yum) 




(cm molec" 1 ) 




2342 


4.27 


C0 2 (v 3 ) 


7.6 xlO" 17 


[1] 


3250 


3.07 


H 2 (n) 


2.0 xlO" 16 


LU 


-2850 


-3.5 


H 2 2 (v 2 + v 6 ) 


5.7 xlO" 17 


[3] 


2139 


4.67 


CO (vi) 


1.1 xlO" 17 


[2] 


2044 


4.89 


CO a (v 3 ) 


8.9 xlO" 17 


[4] 


1307 


7.65 


H2CO3 (C-OHbend) 


1.0 xlO" 16 


[5] 


~ 1040 


-9.6 


3 (v 3 ) 


1.4 xlO" 17 


[6] 



[1] Gerakines et al. 1995; [2] Gibb et al. 2000; [3] Loeffler et al. 
2006; [4] Bennett et al. 2004 [5] Gerakines et al. 2000; [6] Smith et 
al. 1985. 



3. Results 

FigureQJ presents the infrared spectra of H 2 0:C0 2 ice (1:1) at 
13 K, before (top trace) and after different irradiation fluences. 
Each spectrum has been shifted for clarity. The narrow peak at 
2342 cm" 1 is the C0 2 stretching mode (V3). The broad struc- 
tures around 3300 and 800 cm" 1 correspond to the vibration 
modes of water, v\ and vl, respectively. The band at 1600 cm" 1 
is due to the water v 2 vibration mode. The presence of the OH 
dangling band at about 3650 cm" 1 is also observed in the figure 
indicating a high porosity (Palumbo et al. 2006). The inset fig- 
ure shows the newly formed H 2 2 species due to the radiolysis 
of water inside the ice. 

The region between 2200 to 1000 cm" 1 is shown in detail 
in Fig.QJ. The ion fluence of each spectra is also given. In this 
spectral region several peaks grow as a function of fluence: they 
correspond to the appearance of new species formed by the 52 
MeV Ni atoms bombardment, including CO, CO3, O3, H 2 C03, 
H 2 CO (tentative), HCOOH (tentative). Some IR features orig- 
inated from unidentified organic species appear around 1350 
and 1450 cm" 1 . 

The evolution of the column density as function of fluence 
is shown in Fig.flJ. The infrared band strengths used for the de- 
termination of the column density are listed in TableQ] The col- 
umn density of water presents constant concentration of about 
5 x 10 17 molecules cm" 2 , decreasing very slowly as the fluence 
increases. This fact is attributed to an approximate compen- 
sation between persisted water condensation (layering) and its 
disappearance by sputtering and dissociation. 
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H,0:CO, 13 K ice (1:1) 




H,0:CO, 13 K ice (10:1) 



Fig.l. a) Infrared spectra of H 2 0:C0 2 13 K ice (1:1) be- 
fore (top dark line) and after different irradiation fluences. b) 
Expanded view from 2200 to 1000 cm -1 , c) Molecular column 
densities derived from the infrared spectra during the experi- 
ment. The lines are to guide the eyes. 




Fluence (10 ions cm") 



Fig. 2. a) Infrared spectra of H 2 0:C0 2 13 K ice (10:1) be- 
fore (top dark line) and after different irradiation fluences. b) 
Expanded view from 2200 to 1000 cm -1 , c) Molecular column 
densities derived from the infrared spectra during the experi- 
ment. The lines are to guide the eyes. 



The initial enhancement in the water column density maybe 
associated with the compaction effect produced by ion bom- 
bardment as discussed elsewhere (Palumbo 2006, Pilling et al. 
2010). This compaction changes the band strength of some vi- 
brational modes. Leto and Baratta (2003) showed that the band 
strength of water v\ vibrational mode undergoes a 10% increase 
during the first ion impacts (low fluence) on ice in experiments 
employing ions with energy of dozens of keV. Another expla- 
nation may be a strong water layering just on the begging of 
irradiation. The C0 2 abundance reaches its half value at a flu- 
ence of about 4 x 10 12 ions cm 2 . The CO abundance increases 
very rapidly, reaching a maximum around 3-4xl0 12 ions cm" 2 
and decreasing after that. The same behavior was observed 
for the O3, H 2 CC<3 and CO3, all produced by the radiolysis 
of C0 2 inside the ice. For O3, H 2 CC<3 and CO3, the column 



density maximum occurs at 4 xlO 12 , 3 xlO 12 and 1 xlO 12 ions 
cm" 2 , respectively. In the case of H 2 2 the maximum occurs 
at around 1 xlO 12 ions cm" 2 and remains constant (~ 1 xlO 16 
molec cm" 2 ) until the end of irradiation (up to a fluence of 1 
xlO 13 ions cm" 2 ). 

The infrared spectra of H 2 0:C0 2 (10:1) ice before and af- 
ter different irradiation fluences up to 5xl0 12 ions cm" 2 are 
given in Fig. [2^. In contrast to H 2 0:C0 2 (1:1), this ice does not 
present the OH dangling band at about 3650 cm" 1 , may be due 
to the lower deposition rate or due to the low amount of C0 2 
on ice. As suggested by Bouwman et al. (2007), the presence of 
impurities like CO and C0 2 on amorphous water ices increases 
the ice's porosity and the presence of the OH db become more 
prominent. 
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Fig. 3. a) Infrared spectra of pure water ice at 13 K obtained for Fig. 4. a) Infrared spectra of pure C0 2 ice at 13 K after differ- 
different 52 MeV Ni fluences. Inset figures present details of ent 52 MeV Ni fluences. Inset figures present details of newly 
newly formed species, b) Molecular column densities derived formed species, b) Molecular column densities derived from 



from the infrared spectra during the experiment. 



the infrared spectra during the experiment. 



Figure [2J3 shows details of the infrared spectra for regions 
between 2200 to 1000 cirT 1 . The wavenumber of some uniden- 
tified IR frequencies are given. The infrared peaks tentatively 
attributed to H 2 CO (~ 1470 cm 1 ) and HCOOH (~ 1220 cm 1 ) 
seem to be smaller than in the case of H20:C02 (10:1) ice. As 
in the previous case the appearance of the CO v\ line at 2139 
cm -1 is due to the CO2 radiolysis. 

The variations of the column densities of the most abun- 
dant molecules observed during the irradiation of H20:CC<2 ice 
(10:1) by 52 MeV Ni ions as function of fluence are given in 
Fig. [2}:. Peaks due to O3, H2CO3 or CO3 molecules were not 
observed. 

The infrared spectra of pure H2O and CO2 ices at 13 K for 
different 52 MeV Ni fluences are shown in Fig.fJk and Fig.|4^, 
respectively. Upper curves indicate the virgin ice and figure in- 
sets display peak details of the newly formed species. In the 
case of the radiolysis of pure water, only the formation of H2O2 
(vibration mode V2+V6) is observed. In contrast, after the ra- 
diolysis of CO2 pure ice, at least three different compounds 
are formed: CO (vi), O3 (V3) and CO3 (vi) (see also Seperuelo 
Duarte et al. 2009). 

Figure [3J) shows the evolution of water molecular column 
density derived from different vibration modes (vi, vi and v£) 
and the column density of H2O2 as a function of ion fluence. 
The band strengths employed for water vibration modes vi 
(1650 cirT 1 ) and Vl (~ 800 crrr 1 ) were obtained by Gerakines 
et al (1995) and Gibb et al. (2000), respectively. The evolution 



of CO2 column density, derived from different vibration modes 
(V2, V3, vi + V3 and 2V2+V3), is given in Fig.[4}5: the distinct result 
are practically the same, as expected. The band strengths em- 
ployed for CO2 vibration modes V2 (670 cm -1 ), 2v2 + V3 (3599 
cm and vi +v 3 (3708 cirT 1 ), were obtained by Gerakines et 
al. (1995). The column densities of the new species produced 
by radiolysis are also shown and will be discussed further. 

The infrared absorption profile of water (vi) and carbon 
dioxide (V3) in pure and mixed (H20:C02 (1:1)) ices as a func- 
tion of ion fluence are shown in Fig. [5] Upper panels (a and 
c) indicate pure ices. In the case of mixed ice, the selected IR 
features are illustrated in the bottom panels (b and d). In all 
panels the upper curves indicate the non irradiated spectrum. 
Each spectrum has a small offset for clearer visualization. The 
bombardment by heavy ions slightly shifts the water vi vibra- 
tion mode to low frequencies. Since different water clusters are 
responsible for this IR band (Paul et al. 1997) this suggests that 
small water clusters within the ice are being converted to larger 
cluster structures. Also large clusters may be less radiation sen- 
sitive than small clusters. This behavior is enhanced in the case 
of mixed ices (Fig.[5j3). 

Figure [5J5 presents the OH dangling bands (-3670 cm -1 ) 
attributed to water molecules weakly adsorbed into micropores 
inside the ice (see Palumbo 2006). As the fluence increases, 
the micropores collapse and the OH dangling vibration modes 
decrease. We have shown previously (Pilling et al. 2010) that 
the ice compaction produced by heavy cosmic rays are at least 
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3800 3400 3000 2600 2400 2360 2320 2280 2240 



Fig. 5. Selected profile of water and CO2 vibration mode at dif- 
ference fiuences up to 10 13 ions cirT 2 . Pure water v\ vibration 
mode (a), water v\ mode in mixed ices (b), pure CO2 V3 mode 
(c) and CO2 V3 mode in mixed ices (d). See details in text. 

3 orders of magnitude higher than that produced by (0.8 MeV) 
protons. 

The profile of the V3 vibration mode of carbon dioxide 
appears to be very sensitive to the presence of water. It be- 
comes very broad and also presents a shift toward the lower 
frequencies when water is mixed with CO2 (Fig. |5ji). As the 
ion fluence increases, the V3 band of CO2 becomes sharp and 
the peak moves toward higher frequencies becoming similar to 
pure CO2 ice. A detailed discussion on the effects of CO2 and 
H2O on band profiles observed in mixed H20:C02 ices is given 
elsewhere (Oberg et al. 2007). 

4. Discussion 

4.1. Dissociation cross section 

As discussed previously by Seperuelo Duarte et al. (2009) and 
Pilling et al. (2010), the complex interaction between an en- 
ergetic heavy ion with an ice target may be described through 
partial scenarios by considering successive aspects of the phe- 
nomenon. Since the impact of 52 MeV Ni atoms lies in elec- 
tronic energy loss domain (projectile velocity > 0.01 cm //s -1 ) 
most of the deposited energy leads to excitation/ionization of 
target electrons (electronic energy loss regime). In turn, the 
electrons liberated from the ion track core (0.3 nm of diame- 
ter for 52 MeV Ni atoms; Iza et al. 2006), transfer their en- 
ergy to the surrounding condensed molecules (~ 3 nm). The 
re-neutralization of the track proceeds concomitantly with the 
local temperature rise, leading to a possible sublimation. 

As discussed previously (Pilling et al 2010), the column 
density rate (a function of fluence) for each condensed molec- 
ular species during radiolysis can be written as: 

= X °"/.y^ + Li ~ ad ' iNi ~ FA(F) (2) 



where Yjj °~fjjNj represents the total molecular production rate 
of the i species directly from the j species, L, is the layering, 
(Tdj is the dissociation cross section, Y-, is the sputtering yield 
of and £lt(F) is the relative area occupied by the i species on 
the ice surface. 

Considering that the analyzed compounds cannot react in 
one-step process to form another species originally present in 
ice (i.e., 07 ^ ~ 0), in the case of pure ices Eq.[2]can be written 
by: 

Ni = (N ,i - iVoo.j) exp(-cr d jF) + N^j (3) 

where Ncaj = (Li - Yi)jad,i is the asymptotic value of column 
density at higher fiuences due to the presence of layering. Ni 
and No,i are the column density of species i at a given fluence 
and at the beginning of experiment, respectively. In the case of 
pure CO2 the layering yield is negligible Lqoi = 0. 

In the case of water-mixed ices, after the rapid ice com- 
paction phase at the beginning of the irradiation, water layer- 
ing tends to recover the ice surface, Q;(F) — » 6 a , preventing 
progressively sputtering of CO2 molecules. At the end of these 
processes, the system of Eq.|2]is decoupled and can be solved 
analytically by: 

N = (Nq - N x ) exp(-cr^) + N„ for H 2 (4) 
and 

N = No sxp(-cr d F) for C0 2 (5) 

where Noo = (L - Y) / cr d is the asymptotic value of column 
density of water due to the presence of layering. 

As discussed by Loeffler et al. (2005), the radiochemical 
formation yield (G/) of a given compound per 100 eV of de- 
posited energy, at normal incidence, can be expressed in terms 
of the formation cross section (erf) and the stopping power (S), 
in units of eV cm 2 /molec, as: 

Gf = 100— molecule per 100 eV (6) 
S 

This definition can be extended for radiochemical dissoci- 
ation (destruction) yield of a given compound per 100 eV of 
deposited energy (G d ) by replacing the formation cross section 
in Eq.[6]by the negative value of the dissociation cross section 
(-<Td). Therefore, negative Gd values indicate that molecules are 
being dissociated or destroyed after energy deposition into the 
ice. 

Adopting the S values from the Stopping and Ranges mod- 
ule of the SRIM2003 package (Ziegler 2003), the values of the 
radiation yield G in the experiments can be determined and 
compared with literature values. Table 2 presents the radio- 
chemical destruction yields of H2O and CO2 per 100 eV of 
deposited energy in each calculated model. The stopping pow- 
ers of 52 MeV Ni ions for pure water ice and for CO2 ice are 
S = 1.4 X 10~ 12 eV cm 2 /H 2 and 2.6 xl0~ 12 eV cm 2 /C0 2 , 
respectively. In the case of mixed H20:CC»2 ices the stopping 
power is determined by interpolation between the pure ice val- 
ues. We considered that the chemical changes during the irra- 
diation does not affect the stopping power. 
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Table 2. Dissociation cross section, radiochemical yield and sputtering values of H2O and CO2 obtained from the radiolysis of 
pure water ice and mixed \vater-CO2 ices by 52 MeV Ni ions. 



Species 


Mixture 


cr d 


-G rf 


/V„ 


Y 


L* 


N 


Model 




(H 2 0:C0 2 ) 


(IQ- 13 cm 2 ) 


molec /100 eV 


(10 17 molec cm" 2 ) 


(10 4 molec ion" 


) (10 4 molec ion" 


') (10 17 molec cm" 2 ) 




H 2 


(1:0) 


1.1 


7.48 


12 


1° 


14 


15 


1 




(10:1) 


~1 


4.90 


~4 


1" 


~5 


6.9 


2 




(1:1) 


1.0 


6.35 


3.0 


1" 


4.0 


5.9 


3 


C0 2 


(0:1) 


1.8 


6.90 


NA + 


2.2 





9.1 


4 




(1:1) 


1.6 


7.84 


NA 


1.2* 





6.2 


5 a 




(1:1) 


2.1 


10.3 


NA 


C 





6.2 


5b 




(10:1) 


~1 


6.35 


NA 


~0.2 rf 





-0.8 


6a 




(10:1) 


-0.7 


4.45 


NA 


C 





-0.8 


6b 



8 Taken from Brown et al. (1984). b Estimated to be a half of the value determined from pure CO2 ice. c No CO2 sputtering. Assuming that the 
water layering is high enough to fully cover the C0 2 molecules on the surface. d Estimated to be a tenth of the value determined from pure C0 2 
ice. * NA=Non applied. * For H 2 0: Derived from L = N^aj + Y. For C0 2 : Assuming no layering due to residual C0 2 . 



Figure [6] presents the fitting curves for the H2O and CO2 
column densities employing Eq. [3] (for pure ices), and Eqs. [4] 
and [5] (for mixed ices). Numeric labels indicate parameters 
listed in Table[2] The H2O column density data level off around 
N x - 1.2xl0 18 ,4xl0 17 , and 3 xlO 17 molec cm" 2 for pure wa- 
ter, H 2 0:C0 2 (10:1) and H 2 0:C0 2 (1:1) ice samples, respec- 
tively. The average value for the dissociation cross section of 
water, employing Eq.[3] is crv ~ 1 x 10" 13 cm 2 . The sputtering 
yield for water measured by Brown et al. (1984) was extrapo- 
lated for the 52 MeV Ni ions impact as discussed in Pilling et 
al. (2010). The obtained value is Y = 1 x 10 4 molecules per im- 
pact. The estimated average water layering of both experiments 
was within the 4-14xl0 4 molec ion" 1 range, obtained from the 
relation L = Ncacr^ + Y. Due to the large value of water layering, 
CO2 species in the mixed ices were recovered by a H2O film 
and the CO2 ice sputtering yield is considered negligible in the 
current experiments; therefore, the data are adjusted directly by 
Eq.HJ 

The fitting parameters (dissociation cross section, radio- 
chemical yield, sputtering, layering and the initial relative 
molecular abundance of each species in the ices) are listed in 
Table [2] Pure H2O and CO2 ice were fitted by the parameters 
of models 1 and 4, respectively. For water species the layer- 
ing yield was determined by L = N inf +<x rf + Y, For CO2 
species (model 4 to 6b), it is assumed no layering due to resid- 
ual CO2. Model 5a is the fitting for the CO2 species in the 
mixed H20:C02 ice (1:1) and the sputtered value was consid- 
ered to be half the value from pure CO2 ice. 

In model 5b it is assumed that the water layering is high 
enough to fully cover the CO2 on the surface producing a negli- 
gible sputtering yield of CO2 molecules. From the comparison 
of models 5a and 5b we observe that the error in the dissoci- 
ation cross section is lower than 30%. Model 6a concerns the 
fitting for the CO2 species in the mixed H20:C02 (10:1) ice; 
the sputtering yield is considered to be a tenth of the value de- 
termined from pure CO2 ice. 

From models 5 we observe that the precise determination of 
sputtering is possible only for the ices monitored at higher ion 
fluences (~ 10 13 ion cm" 2 ) since different sets between sput- 
tering and dissociation cross section can be adjusted for the 



data with low fluences (< 5 x 10 12 ions cm" 2 ). Due to the de- 
generacy of such model for low ion fluence IR data, the error 
observed in the dissociation cross section on models 6 is about 
40%. Table [2] shows that the dissociation cross section at typi- 
cal astrophysical ice ([CC^l/tFhO] - 0.05-0.1) are similar and 
equal to - 1 x 10" 13 cm" 2 . 
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Fig. 6. Variation of the experimental column densities of H2O 
and CO2 as a function of fluence. Lines represent the fittings 
using Eq.[3]for pure ices and Eqs. [4] and [5] for mixed ices. The 
parameters models are listed in Table [2] 



4.2. Synthesis of new species 

The evolution of the column density of the newly formed 
species from H20:C02 ices as a function of fluences is shown 
in Fig. [7^ and Fig.lTJ). The ratio of H2O2 column density over 
its parent molecule (H2O) initial column density, Nh-,0-, /Wo,h 2 o> 
as a function of fluence in 3 different water-concentration ices 
(pure H2O ice, -50% and -90%) is given in Fig. [7^. The 
evolution of newly formed species CO, O3, and H2CO3 from 
CCh-rich ices as a function of fluence is presented in Fig. [7}}. 
Nt/No,co 2 indicates the column density ratio of a given pro- 
duced species over its parent molecule (CO2) initial column 
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Table 3. Formation and dissociation cross section of newly formed species from the radiolysis of H20:C02 ices by 52 MeV Ni 
ions. 



Species 


Parental 
species 


Mixture 
[H 2 0:C0 2 ] 


°7 

(10- 13 cm 2 ) 


(10- 13 cm 2 ) 


Gf 

(molec /100 eV) 


-G d 
(molec /100 eV) 


Model 


H 2 2 


H 2 


[1:0] 


0.19 


10 


1.3 


68.0 


1 




H 2 


[10:1] 


0.13 


9.3 


0.87 


59.1 


2 




H 2 


[1:1] 


0.10 


6.9 


0.48 


33.8 


3 


CO 


C0 2 


[0:1] 


1.8 


11 


6.8 


42.1 


4 




C0 2 


[1:1] 


1.2 


7.3 


5.8 


35.8 


5 




C0 2 


[10:1] 


0.91 


4.4 


5.8 


27.9 


6 


O, 


C0 2 


[0:1] 


0.27 


16 


1.1 


61.3 


7 




C0 2 


[1:1] 


0.06 


5.0 


0.30 


24.5 


8 


co 3 


C0 2 


[0:1] 


0.09 


97 


0.33 


372 


9 




C0 2 


[1:1] 


0.01 


31 


0.06 


152 


10 


H2C03 


H 2 and C0 2 


[1:1] 


0.03 


9.9 


0.14 


48.6 


11 



density. Three different CC^-rich ices are analyzed (pure CO2 
ice, ~50% and ~10%). The data was fitted by the simple expo- 
nential associate expression: 



Ni/N , p = At + Bi[l - expC-o-^F)] 



(7) 



where Ni and No, p represents the column density of a given 
produced species i over its parent molecule p (H2O or CO2). 
A,- and B, are constants which indicate the initial and the max- 
imum amount (after radiolysis) of species i on the ice, oy ,■ is 
the formation cross section and F is the ion fluence. 

As the fluence increases, the number of produced 
molecules also increases and reaches a maximum. At this 
moment, the number of molecules produced by radiolysis of 
parental species is equal to the number of new molecules that 
dissociate possibly in the reverse reaction set (e.g. p ^± i + j + 
k+...). At the equilibrium the column density of the parental 
and daughter species, N eq , p and N eq j, respectively, are related 
by the expression: 



Neq,p O-fj = N e qJ CT d j 



(8) 



where 07 , and o~d,i is the formation cross section and dissocia- 
tion cross section of daughter species i. 

Since the column densities of daughters species i are much 
smaller than the one observed for parental species p, as a first 
approximation we have the following relation between them: 



N„, 



N,„ 



1 



(9) 



Noj, N , p 

During the whole irradiation time, when a new species is 
being formed from the radiolysis of its parental species, a frac- 
tion of them is also being dissociated and possibly being sput- 
tered from the ice. This last effect is not being considered in this 
model. Since not all the parental molecules are converted into a 
given species, some could be sputtered from the ice and others 
can react to produce different daughter species, the amount of 
parental molecule always decreases. Consequently, the column 
density of a given daughter species decreases after it reaches a 
maximum (equilibrium stage). This behavior can be observed 
in the case of the evolution of CO3 species as a function of 
fluence (Fig.|7p). 



Calling the maximum amount of newly formed species i at 
the equilibrium as Nooj and using Eq|9l we can be express the 
constant B, of Eq. |7]by: 

(10, 



d; = tt — - I - 



10 



HP, (from H,0 pure) « 




H.O n (from H.O at -90%) 
H,0, (from H 2 at 50%) 



a) 



» t 

CO (from C0 2 at 50%) 

O, (from CO pure) 



O. (from CO, at 50%) 

-3 a 

H 2 CO 3 (fromCO,at50%) 

, CO (from CO, pure) 



b) 



2 3 
Fluence (10 l_ ions cm" ! ) 



Fig. 7. The evolution of newly formed species from the bom- 
bardments of H20:CC»2 ices by 52 MeV Ni as a function of 
fluence. a) Production of H2O2 from H2O molecules in dif- 
ferent ices, b) Production of CO, O3, CO3 and H2CO3 from 
CO2 molecules in different ices. The relative abundance of each 
parental species are indicated. The lines indicate the fittings us- 
ing Eq.QT]and the model parameters are given in Table [3] 
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Moreover, considering that the number of initial daughter 
species i is negligible (A, = 0) and, after the combination of 
Eq.|9]with Eq.[l0] we rewrite Eq.|7]as: 

where (Tdj and cr/j are the dissociation and formation cross 
section of species i. 

The lines observed in the Fig. [7^ and Fig. |7J) were ob- 
tained by fitting the experimental data with EqQT| The pro- 
posed model is valid from the beginning of irradiation up to 
the fiuence in which the column density reaches its maximum 
value. This indicates the equilibrium point between formation 
and destruction promoted by radiolysis. The formation and dis- 
sociation cross section as well respective radiolysis yield deter- 
mined for the newly formed species are given in Table [3] We 
observe that CO presents the highest formation cross section 
among the studied newly formed species. 

4.2.1. Hydrogen peroxide, H 2 2 

In recent years, different groups have studied the formation of 
hydrogen peroxide by ion bombardment of pure water ice and 
mixed water-CC<2 ices (Moore & Hudson 2000, Strazzulla et 
al. 2003, 2005a, 2005b; Baragiola et al. 2004; Loeffler et al., 
2006; Gomis et al. 2004a, 2004b). H 2 2 was also observed 
from fast electron bombardment of pure water ices (Baragiola, 
et al. 2005; Zheng et al. 2006) and also from low energy (3-19 
eV) electron bombardment (Pan et al. 2004). 

Following Teolis et al. (2009), the most commonly H2O2 
formation mechanisms involves the reaction between two OH 
radicals coming from the radiolysis of water molecules, as 
given by: 

H 2 ^ H 2 0* — >OH + H (12) 
OH + OH^H 2 2 (13) 

where CR denotes cosmic rays. 

These OH radicals are strongly hydrogen bonded to water 
molecules (Cooper et al., 2003), and it is not until above 80 
K that they can diffuse within a water-ice lattice (Johnson and 
Quickenden, 1997). As pointed out by Cooper et al. (2008) en- 
ergetic OH radicals can diffuse short distances along ion tracks 
and react at 10 K, but bulk diffusion probably does not occur. 

The decrease in the H2O2 production with temperature in- 
crease was investigated by Zheng et al. (2006) and Loeffler et 
al. (2006). This indicates that electron scavenging may play a 
critical role in the radiation stability of H2O2 in pure water- 
ice experiments. Moreover, Gomis et al. (2004a) found that the 
H2O2 yields were dependent on projectile: ices irradiated with 
low energy (30 keV) C + , H + and + ions produced more H2O2 
at 77 K than 16 K, while N + and Ar + had no temperature de- 
pendence on the H2O2 yield. The authors also observe that the 
H irradiation produces a much lower quantity of H2O2 than for 
the other heavy ions suggesting that the energy deposited by 
elastic collisions plays an important role in such a process. 



As discussed by Zheng et al. (2006) and Loeffler et al. 
(2006), there is another reaction sequence to produce hydro- 
gen peroxide in astrophysical ices from the addition of oxy- 
gen atom (a product of the by the dissociation of oxidant com- 
pounds such as O2, CO, CO2, H2O, etc.) to water molecule: 

O + H 2 H 2 0-0 'Tfff " h 2 + + O " -> H 2 2 (14) 

This reaction route involves an extra ionization stage triggered 
by CR, UV/X-ray photons or fast electrons to produce H2O2 
(Loeffler et al. 2006), which implies a quadratic dependence on 
irradiation fiuence being a secondary order reaction processes. 
In the present work we consider negligible the H2O2 formation 
via oxygen addition to H2O. Future investigations with isotopic 
labeling could help to clarify this issue and quantify which frac- 
tion of O atoms produced from the dissociation of CO2 or H2O 
may react with H2O to form H2O2. 

Table [4] presents the formation and dissociation cross sec- 
tion, as well the radiochemical yield, of H2O2 obtained from 
the processing of pure water ice and mixed H20:C02 ice. Both 
the formation and dissociation cross section of H2O2 decrease 
as the relative abundance of H2O in the ice decreases, a con- 
sequence of larger averaged distance between parental species. 
From Fig. [7] we also observe that the presence of CO2 inside 
the ice decreases the H2O2/H2O ratio. This points in the op- 
posite direction than the results observed from the irradiation 
of water-C02 by light ions (Strazzulla et al. 2005b; Moore and 
Hudson 2000). From Table |4] we observe that the formation 
cross section of H2O2 from the bombardment of pure water ice 
by heavy and energetic ions is 10 5 higher than the value ob- 
tained by the impact of light or/and slower projectiles. 

The asymptotic value for the H2O2/H2O ratio presents a 
slight reduction with the enhancement of initial CO2 in the ice. 
This clearly shows that H2O2 formation via O + H2O in which 
the oxygen comes from another oxidant compound, such as 
CO2, is indeed a secondary-order process (Loeffler et al. 2006). 
Our value is about 8 times higher than in the case of irradiation 
whit 30 keV protons, but very similar to the asymptotic values 
obtained after the impact with other 30 keV ions such as C + , N + 
and + (Gomis et al 2004a). A comparison between our radio- 
chemical yield and cross sections with literature have revealed 
higher values for both the formation and destruction of H2O2 
during bombardments of ices by heavy and energetic ions. In 
other words the chemical processing is very enhanced by the 
processing of ices by heavy and energetic and highly charged 
ions in comparison with light ions. 

Hydrogen peroxide has been found on the surface of 
Europa by identifying both an absorption feature at 3.5 /mi 
in the Galileo NIMS spectra and looking at the UV spectrum 
taken by the Galileo ultraviolet spectrometer (UVS) (Carlson 
et al. 1999). Following the authors, the relative abundance with 
respect to water in Europa is H 2 2 /H 2 ~ 0.001 (=0.1%). As 
the Jupiter satellite is immersed into the intense magnetosphere 
of the planet, the authors also suggest that radiolysis is the 
dominant formation mechanism of such a molecule. As dis- 
cussed by Hendrix et al. (1999) this compound has also been 
detected on the other icy Galilean moons of Jupiter in the ul- 
traviolet, such as Callisto and Ganymede, although it is less 
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Table 4. Comparison between the formation and dissociation cross section of H2O2 from the processing of pure H2O ices and 
H20:C02 ices. The formation and dissociation radiochemical yield of H2O2 as well the asymptotic H2O2/H2O ratio in each 
experiment are also given. 



Ices 


Temp. 


Projectile 


°7 




G f 


-G rf 


I^Oj/I^O 1 


Ref. 




(K) 


(Energy) 


(10" 13 cm 2 ) 


(10" 13 cm 2 ) 


(molec / 100 eV) 


(molec / 100 eV) 


(%) 




H 2 


13 


Ni 13+ (52MeV) 


0.19 


10 


1.3 


68.0 


1.9 


[1] 


H 2 0:C0 2 (10:1) 


13 


Ni 13+ (52MeV) 


0.13 


9.3 


0.87 


59.1 


1.6 


[1] 


H 2 0:C0 2 (1:1) 


13 


Ni 13+ (52MeV) 


0.10 


6.9 


0.48 


33.8 


1.5 


[1] 


H,0 


16 


H + (30 keV) 


~ 0.001 


0.06 


0.45 


271 


0.21 


[2] 


H 2 


16 


C + (30 keV) 


~ 0.001 


0.07 


0.23 


15.8 


1.7 


[2] 


H 2 


16 


N + (30 keV) 


~ 0.001 


0.09 


0.19 


17.1 


1.6 


[2] 


H 2 


16 


0+ (30 keV) 


0.002 


0.13 


0.34 


22.3 


1.4 


[2] 


H 2 


16 


Ar + (30 keV) 


0.002 


0.07 


1.34 


46.8 


2.2 


[2] 


H,0 


20 


H + (100 keV) 


~ 0.001 


0.20 


0.38 


76.1 


0.75 


[3] 


H 2 


16 


H + (200 keV) 


~ 0.001 


0.09 


0.49 


43.9 


1.2 


[4] 


H 2 


16 


He + (200 keV) 


0.003 


0.13 


0.37 


16.1 


2.4 


[4] 


H 2 


16 


Ar 2+ (400 keV) 


-0.001 


0.04 


0.07 


2.79 


6.0 


[4] 



Asymptotic value; [1] this study; [2] Gomis et al. 2004a; [3] Loeffler et al. 2006; [4] Gomis et al. 2004b. 



evident on the darker satellites than on Europa. The estimated 
ratio on Ganymede is H 2 2 /H 2 ~ 0.003 (=0.3%) (Hendrix et 
al. 1999). 

H2O2 has also been suggested to be present in icy mantles 
on grains in dense clouds as pointed out Boudin et al. (1998). 
Based on observations and laboratory studies, they have esti- 
mated the relative abundance with respect to water at the dense 
molecular cloud NGC 7538:IRS 9 of H 2 2 /H 2 ~ 0.05 (=5%). 
Following the authors, hydrogen peroxide can be also produced 
on grain surfaces via the hydrogenation of molecular oxygen. 

4.2.2. Ozone, 3 

Ozone has been detected after irradiation of mixed H20:C02 
(1:1) ices and pure CO2 ice at 16 K by 1.5-200 keV ions (pro- 
tons and He + ) (Strazzulla et al. 2005a, 2005b). According to 
these authors, the O3 production from radiolysis is highly atten- 
uated in the case of CC>2-poor mixed ices, e.g. for H2O/CO > 
2.5. We do not observe ozone among the compounds produced 
by radiolysis of H2:C02 (10: 1) ice. The presence of ozone has 
also been observed after the irradiation of CO:C<2 (1:1) ice by 
60 keV Ar ++ and 3 keV He + (Strazzulla et al 1997). The au- 
thors suggest that the detection of ozone and the absence of 
suboxides (e.g. C2O and C3O) in interstellar grains would indi- 
cate a dominance of molecular oxygen in grain mantles. These 
suboxides have not been detected in our experiments. 

Seperuelo Duarte et al. (2010) have bombarded pure CO 
ice with 50 MeV Ni +13 and observe ozone and several subox- 
ides (C2O, C3O2, C5O2) among the radiolysis products. They 
determined the ozone formation cross section is 3 xlO" 16 cm 2 . 
Ozone and suboxides were also observed after irradiation of 
pure CO ice at 16 K by 200 keV protons (Palumbo et al. 2008). 

We have detected ozone after the radiolysis of pure CO2 
ice and mixed H 2 0:C0 2 (1:1) ice by 52 MeV Ni +B . The deter- 
mined formation cross section of ozone is cr 03 = 0.27 x 10" 13 
cm 2 and 0.06 x 10~ 13 cm 2 in pure CO2 ice and H 2 0:C02 ice, 
respectively. These values indicate that the O3 formation cross 



section decreases as the relative abundance of the parental CO2 
compound in the ice declines. The value obtained for ozone 
formation via pure CO2 is 2 times higher than the value deter- 
mined previously in similar radiolysis experiments involving 
pure 18 C0 2 irradiated by 46 MeV 58 Ni +n (Seperuelo Duarte 
et al. 2009). 

In the solar system, ozone has been observed on Ganymede 
(Noll et al. 1996) and on Saturn's satellites Rhea and Dione 
(Noll et al. 1997). Following Strazzulla et al. (2005a) ozone 
could be formed where fresh CO2 rich layers are exposed to ra- 
diation. Due to the strong silicate absorption band around 1040 
cm" 1 the observation of ozone in interstellar medium is very 
difficult. However, the appearance of ozone seems conceivable, 
since it is observed in laboratory experiments involving the pro- 
cessing of astrophysical ices analogs. 

4.2.3. C0 3 

This compound has been observed after the UV processing 
of pure C0 2 , CO and 2 10 K ices (Gerakines et al. 1996; 
Gerakines & Moore 2001) and also on radiolysis of pure CO2 
by 0.8 MeV protons (Gerakines & Moore 2001) and 1.5 keV 
protons (Brucato et al. 1997). The bombardment of mixed 
H 2 0:C0 2 ices by 3 keV He + , 1.5 keV and 0.8 MeV protons 
also shows the formation of CO3. Furthermore, H2CO3 and 
possibly O3 are produced (Brucato et al. 1997; Moore and 
Khanna 1991). 

Seperuelo Duarte et al. (2009) have observed this species 
in the radiolysis of of C 18 2 by 46 MeV 58 Ni 11+ up to a fi- 
nal fluence of 1.5 x 10 13 cm" 2 . In our measurement the forma- 
tion cross section of CO3 from radiolysis of pure CO2 ice is 
0.08 x 10" 13 cm" 2 , roughly 5 times lower than the value ob- 
tained previously by Seperuelo Duarte et al. (2009). This value 
presents a decrease in a presence of water in the ice, as ob- 
served in the case of H20:C02 (1:1) ice (see Table[3]). The for- 
mation cross section of CO in our experiments on pure CO2 is 
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roughly the same than the one determined by Seperuelo Duarte 
et al. (2009). 

The CO3 destruction cross section due to radiolysis is very 
large when compared to the other compounds formed into the 
ices. To increase the accuracy of the cross sections next mea- 
surements should cover large data sets inside the fiuence range 
between 1 - 1 Ox 1 1 1 ions cm" 2 . On the contrary of UV photoly- 
sis of pure CO2 ices (Gerakines et al. 1996) we do not observed 
the formation of C3O at 2243 cm" 1 from the radiolysis of CO2 
ices. 

Despite of the prediction of CO3 among the compounds 
in the processed astrophysical ices (e.g. Elsila et al. 1997, 
Allamandola et al. 1999), this species has not been detected 
conclusively in space yet (Elsila et al. 1997, Ferrante et al. 
2008). 

4.2.4. H2CO3 

Carbonic acid (H2CO3) has been observed in several experi- 
ments involving ion bombardment of H20:CC<2 ices (Moore & 
Khanna, 1991; DelloRusso et al. 1993; Brucato et al. 1997). It 
was also observed after the radiolysis of pure CO2 ice at 10 K 
by 1.5 keV protons (Brucatto et al. 1997), indicating that im- 
planted hydrogen ions are incorporated in the target to form 
new bonds to produce H2CO3. 

As discussed by Gerakines et al. (2000) the formation of 
H2CO3 from ion bombardment of H20:C02 ices is ruled by 
two main reaction schemes. First, the direct dissociation of 
H2O molecules by incoming projectile: 

H 2 — > OH + H + + e" (15) 

where products such as H2, H2O2; H30 + , and HO2 are even- 
tually formed by reactions involving the primary dissociation 
products. The next step is the electron attachment on CO2 or 
OH producing reactive compounds which quickly react with 
each other to produce bicarbonate (HCO3 ): 

COj+OH^HCOj and/or C0 2 + OH ^HCO" (16) 

Finally, bicarbonate reacts with a proton to produce H2CO3: 

HCO3 + H + ^H 2 C0 3 (17) 

From the evolution of the IR peak at 1307 cm" 1 during 
the radiolysis of H 2 0:C02 (1:1) ice by 52 MeV Ni ions we 
have determined the formation and the dissociation cross sec- 
tion of H2CO3 and the values obtained are 3.0 x 10" 15 cm 2 
and 9.9 x 10" 13 cm 2 , respectively. This formation cross section 
is 8 orders of magnitude higher than the one derived from the 
UV photolysis of H 2 0:CO (1:1) ices at 18 K (Gerakines et al. 
2000) indicating that heavy ion processing of astrophysical ice 
analogs is very efficient to form H2CO3 compared to UV pho- 
tons. 

H2CO3 is thermally stable at 200 K, higher than the sub- 
limation temperature of H2O. As discussed by Peeters et al. 
(2008), ices containing both H2O and CO2 have been found on 
a variety of surfaces such as those of Europa, Callisto, Iapetus, 
and Mars, where processing by magnetospheric ions or the so- 
lar wind and energetic solar particles may lead to the formation 
ofH 2 C0 3 . 



The astrophysical significance of solid carbonic acid have 
been extensively discussed elsewhere (e.g. Khanna et al. 1994, 
Brucato et al 1997, Strazzullla et al. 1996). H 2 C0 3 has been 
suggested to be present at the surface of several solar system 
bodies with sufficient amounts of CO2 such as Galilean satel- 
lites of Jupiter (Wayne 1995, McCord et al. 1997, Johnson et 
al. 2004) and comets (Hage et al. 1998). Europa, Ganymede, 
and Callisto all exhibit high surface CO2 abundances and these 
satellites are heavily bombarded by energetic magnetospheric 
particles, galactic cosmic rays and solar radiation. 

It has been suggested that the 3.8 fim feature in NIMS spec- 
tra of Ganymede and Callisto arise from H2CO3 (Johnson et 
al. 2004; Hage et al. 1998). Following Carlson et al. (2005) 
the relative abundance of carbonic acid with respect to CO2 in 
Callisto is H2CO3/CO2 ~ 0.01. 

The presence of carbonic acid in interstellar ices seems also 
to be very likely (Whittet et al. 1996) not only in the solid phase 
but also in gas phase (Hage et al. 1998). In the case of the qui- 
escent cloud medium toward the background field star Elias 
16, the expected H2CO3 column density could be as much as 
4.6 x 10 16 molec cm" 2 , corresponding roughly 1% of the CO2 
or 0.2% of the total H 2 in this line of sight (Whittet et al. 
1998). 

As discussed by Zheng & Kaiser (2007), in the history of 
Mars there might be a high concentration of carbonic acid pro- 
duced by radiolysis of surface H20:CC>2 ices. Without a dense 
atmosphere and magnetic field, Mars lacks the power to attenu- 
ate penetrating energetic solar particles, energetic cosmic rays 
or any type of high energy particles. Following Westall et al. 
(1998), if available in sufficient concentrations, carbonic acid 
could potentially dissolve metal ores and catalyze chemical 
reactions and its presence on Mars may also lead to the ex- 
istence of limestone (CaC03), magnesite (MgCC>3), dolomite 
(CaMg(C0 3 ) 2 ) and siderite (FeC0 3 ). 

5. Astrophysical implication 

Figure[8]presents the total flux of heavy ions (12 < Z < 29) with 
energy between 0.1-10 Me V/u inside solar system as a function 
of distance to the Sun. Both Galactic cosmic rays (GCR) and 
energetic solar particles are displayed. The integrated flux of 
heavy and energetic solar ions (black square) at Earth Orbit 
were measured by Mewaldt et al. (2007). Starting from the so- 
lar photosphere abundance and the asumtion that elements with 
first ionization potential of less than 10 eV are more abundant 
in the energetic solar particles than in photosphere by a fac- 
tor of 4.5 (Grevesse et al. 1996) at Earth orbit, the integrated 
flux of heavy and energetic solar particles (4>hsw) with ener- 
gies around 0.1-10 Me V/u was found to be about 1.4 x 10" 2 
cm" 2 s" 1 (black square in Fig. [8). 

The heavy ion integrated flux at Earth orbit from galactic 
sources (heavy cosmic rays) was calculated by using the lunar- 
GCR particle model of Reedy and Arnold (1972) (see also Fig. 
3.20 of Vaniman et al. 1991) and taking into an account the 
relative elemental abundances of 12 < Z < 29 atoms at 1 AU 
(Simpson 1983; Drury, Meyer, Ellison 1999). The value ob- 
tained was ~ 7xl0" 3 cm" 2 s" 1 (blue circle in Fig. [HJ. Assuming 
that the fluxes of solar wind and energetic solar particles are in- 
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Distance to the Sun (AU) 



Fig. 8. Estimated value of the integrated flux of heavy ions 
(12 < Z < 29) with energy between 0.1-10 MeV/u inside so- 
lar system and at interstellar medium as a function of distance 
to the Sun. Both Galactic cosmic rays and energetic solar par- 
ticles are displayed. Square: integrated flux of energetic solar 
particles. Circle and triangle: integrated flux of cosmic rays. 

versely proportional to the squared distance, this value can be 
determined as a function of distance up to the heliopause (~100 
AU) (i.e. the boundary between the solar wind domain and the 
interstellar medium). 

The estimation of the average heavy nuclei cosmic ray flux 
(>Phcr) m interstellar medium was performed by (Pilling et al. 
2010), which considered the value of <Phcr ~ 5 x 10~ 2 cirT 2 s _1 
(blue triangle in Fig.[8]l. This value was rather the same as at the 
outer border of heliopause. An error region (dashed area) was 
introduced in the Fig.[8]to take into an account the uncertainty 
of these estimations. 

The value of integrated flux of heavy ions with energies 
around 0. 1-10 MeV/u from galactic sources and energetic solar 
particles are comparable at Mars orbit (~1.5 AU). However at 
the orbits of Jupiter, Saturn, Uranus, Pluto and at the heliopause 
border, 4>hcrI4>hsw ~ 2 x 10 1 , 1 x 10 2 , 7 x 10 2 , 3 xlO 3 and 2 
xlO 4 , respectively. The total flux of heavy particles with energy 
between 0.1-10 MeV/u at Oort cloud distance was assumed to 
be the same as expected in the interstellar medium. 

Considering the estimated flux of heavy particles of ener- 
getic solar particles {<Phsw) and of heavy nuclei cosmic rays 
(4>gcr) as we H the determined dissociation cross section (ov), 
we can calculate the typical molecular half-lives of frozen 
molecules in astrophysical surfaces due the presence of heavy 
particles by the expression (see Pilling et al. 2010): 



(<f>HSW + $HCR) X <T d 

It is worth to note that the present experiments were per- 
formed at low temperature of 1 3 K, the temperature that is ad- 
equate to ices in the interstellar medium. Therefore it is useful 
to make a comparison between IR spectra of interstellar and 




1800 1700 1600 1500 1400 1300 1000 
Wavenumber (cm 1 ) 

Fig. 9. Comparison between IR spectra of interstellar and lab- 
oratory ices. The top three curves are infrared spectra of young 
stellar sources obtained by the Infrared Space Observatory 
(ISO). Lower traces indicate different laboratory spectra of irra- 
diated H20:C02 ices: a and c (this work); b (Hudson & Moore 
2001); d (Gerakines et al. 2000); e (Strazzulla et al. 2005b); f 
(UV photons, Gerakines et al. 2000). Vertical dashed lines in- 
dicate the location of vibrational modes of frozen H2CO3. 

laboratory ices is shown in Figure [9] The top three curves are 
IR spectra of young stellar sources obtained by the Infrared 
Space Observatory (ISO). The two bold (blue) curves are our 
data. The six bottom curves indicate different laboratory spec- 
tra of processed H20:C02 ices: a) mixture (10:1) irradiated by 
52 MeV Ni ions (this work); b) mixture (10: 1) irradiated by 0.8 
MeV protons (Hudson & Moore 2001); c) mixture (1:1) irradi- 
ated by 52 MeV Ni ions (this work); d) mixture (1:1) irradiated 
by 0.8 MeV protons (Gerakines et al. 2000); e) mixture (1:1) 
irradiated by 0.8 MeV protons (Strazzulla et al. 2005b) and f) 
mixture (1:1) irradiated by UV photons (Gerakines et al. 2000). 
Vertical dashed lines indicate the location of vibrational modes 
of frozen H2CO3 (Moore & Khanna 1991). The presence of 
H2CO3 seems to be independent of the ionization source. 

The three bumps observed on the IR spectra of the young 
stelar objects at about 1600 cm -1 , 1400 cirT 1 and 1350 cirT 1 
(see arrows on Pig[9j» present a good similarity with the features 
present in the IR spectra of H20:CC>2 (10:1) irradiated by 52 
MeV ions (curve a). Nevertheless, up to know, its molecular 
assignment remains unknown. Although very attenuated, these 
features are still observed in the IR spectra of H20:CC>2 (1:1) 
(curve c). This suggests that heavy cosmic rays can be a good 
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candidate to explain some features observed in the IR spectra 
of some interstellar sources. 

The comparison between the current results with the pre- 
vious one (Pilling et al. 2010) reveals that, independently of 
ice constitution (involving H2O, CO2, CO and NH3), the dis- 
sociation cross section due to heavy and energetic cosmic ray 
analogs are in the same range, cr d = 1-2 xl0~ 13 cm 2 . This sup- 
ports the extension of our previous estimative for the half-life 
of ammonia-containing ice, T1/2 = 2-3 x 10 6 years, to all kind 
interstellar grains inside of dense interstellar environments in 
the presence of a constant galactic heavy cosmic ray flux. 

The temperature of solar system ices is about 80 K, about 
5-8 times higher than observed in interstellar ices, therefore 
some molecular species (the most volatile as O3 and CO) are 
not efficiently trapped/adsorbed in/on the ices. This issue can 
makes an enormous difference on surface and bulk chemistry. 
Future experiments employing H20:C02 ices at 80 K will be 
performed to investigate this issue. 

6. Summary and conclusions 

The interaction of heavy, highly charged and energetic ions 
(52 MeV 58 Ni 13+ ) with pure H 2 and C0 2 ices and mixed 
(H20:C02) ices was studied experimentally. The aim was to 
understand the chemical and the physicochemical processes in- 
duced by heavy cosmic rays inside dense and cool astrophysi- 
cal environments such as molecular clouds, protostellar clouds 
as well at the surfaces solar system ices. Our main results and 
conclusions are the following: 

1 . In all experiments containing CO2 (pure ice and mixtures) 
after a fluence of about 2-3 xlO 12 ions crrT 2 , the CO2/CO 
ratio became roughly constant (~ 0.1), independently on 
initial CO2/H2O. The experiments suggest that the abun- 
dances of CO3, O3, H2CO3 in typical astrophysical ices 
([C0 2 ]/[H 2 0] ~ 0.05-0.1) irradiated by heavy ions, should 
be very low, except for those ices with peculiar enrichment 
ofC0 2 . 

2. After a fluence of about 2-3 x 10 12 ions crrT 2 the H 2 2 /H 2 
ratio stabilizes at ~ 0.01, not depending of initial H2O rel- 
ative abundance and column density. 

3. A comparison between our radiochemical yield and cross 
sections with literature have revealed higher values for both 
the formation and destruction of H2O2 during bombard- 
ments of ices by heavy and energetic ions. In other words 
the chemical processing is very enhanced by the processing 
of ices by heavy and energetic and highly charged ions in 
comparison with light ions. 

4. The dissociation cross section of pure H2O and CO2 ices 
are 1.1 and 1.9 xl0~ 13 cm 2 , respectively. In the case of 
mixed H20:C02 (10:1) the dissociation cross section of 
both species is roughly 1 xl0~ 13 cm 2 . 

5. Some IR features observed in ices after the bombardment 
by heavy ions present great similarity with those observed 
at molecular clouds, suggesting that heavy cosmic rays play 
an important role in the processing of frozen compounds in 
interstellar environments. 
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